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ABSTRACT 



Context. A 1000-s flare from a new hard X-ray transient, IGR J08408-4503, was observed by INTEGRAL on May 15, 2006 during 
the real-time routine monitoring of IBIS/ISGRI images performed at the INTEGRAL Science Data Centre. The flare, detected during 
a single one-hour long pointing, peaked at 250 mCrab in the 20-40 keV energy range. 

Aims. Multi-wavelength observations, combining high-energy and optical data, were used to unveil the nature of IGR J08408-4503. 
Methods. A search in all INTEGRAL public data for other bursts from IGR J08408-4503 was performed, and the detailed analysis of 
another major flare is presented. The results of two Swift Target of Opportunity observations are also described. Finally, a study of the 
likely optical counterpart, HD 74 194, is provided. 

Results. IGR .108408^1503 is very likely a supergiant fast X-ray transient (SFXT) system. The system parameters indicate that the 
X-ray flares are probably related to the accretion of wind clumps on a compact object orbiting about 10 13 cm from the supergiant 
HD 74 194. The clump mass loss rate is of the order of 10~ 6 Mq yr~' . 

Conclusions. Hard X-ray flares from SFXTs allow to probe the stellar winds of massive stars, and could possibly be associated with 
wind perturbations due to line-driven instabilities. 

Key words, gamma rays: observations - stars: supergiants - X-rays: individuals: IGR J08408-4503 - stars: early-type - instabilities 



1. Introduction 

A new hard X-ray transient, IGR J08408^4503, was discovered 
by INTEGRAL on May 15, 2006 during the real-time routine 
monitoring of IBIS/ISGRI 20-40 keV images performed at the 
INTEGRAL Science Data Centre (ISDC). It was detected during 
a one-hour long pointing, and further analysis of the lightcurve 
with a 100-s binning showed that the source flared for ~1000 s 
starting at 2006-05-15T18:25:33 with a peak flux o f 250 mCrab 
in the 20^-0 keV energy range ( Gotz et aL||2006| l. The source 
was also detected by the JEM-X instrument, providing a position 
of : RA(J2000) = 08 h 40 m 48 5 .7, Dec(J20 00) = -45°03'41" (+46" 
at 90% confidence) ( |Brandt et al.|2006] >. 

Analysis of the public INTEGRAL archive shows that there 
was at least another outburst from this source on July 1, 2003 
(Mereghetti et al. 2006 ). A detailed analysis of this first flare is 
presented in Sect. |2.1| 

A first Swift target of opportunity observation was carried 
out on IGR J08408-4503 on May 22, 2006 ( [Kennea & Campana 
2006). The XRT observation (in photon counting mode) revealed 
the presence of a faint source in the JEM-X error box at the posi- 
tion : RA(J2000) = 08 h 40 m 48 m , Dec(J2000) = -45°03'30" (with 
an error radius of 5 '.'4 at 90% confidence). Sect. 2.2 gives more 
details about this observation and a following one. 



Preliminary studies of ESO archival optical spectra of the 
star HD 74 194, which is lying at the Swift/XRT position, were 
performed by Mas etti et al.| (|2006), who pointed out that this star 
is likely associated with IGR J08408^503, on the basis of the 
detection of Ha in emission. A deeper analysis of these results 
is presented in Sect. [3] 

Finally, Sect. [4] provides a discussion on the nature of this 
new source, and derives likely system parameters based on the 
observations. 

2. The high-energy source 

2.1. INTEGRAL observations 



All available public INTEGRAL (Win kler et al.|2003| pointings 
less than 10° away from the source position were analysed to 
look for other flares from the source. There are two major peri- 
ods during which the source was frequently observed, summa- 
rized in Table [T] The data were analyzed using OS A^] version 
5.1. A lightcurve at the source position was extracted, with a 
binning of 100 s. 

The detailed study of this lightcurve showed that there is 
only one significant flare, around 2003-07-01T20:03:02. A zoom 
on the ISGRI and JEM-X lightcurves around this flare is shown 
in Fig. [T] together with the hardness ratio. The source was sig- 
nificantly detected with ISGRI for 12.6 ks with a maximum 
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1 The offline scientific analysis (OSA) software is available from the 
ISDC website : http : //i sdc . uni ge . ch 
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Table 1. Summary of the main INTEGRAL public data available 
for IGR J08408-4503. 



Rev. 


Period [MJD] 


Date 


81-88 


52 803-52 827 


12 Jun-6 Jul 2003 


137-141 


52970-52 984 


27 Nov-1 1 Dec 2003 



flux reaching 11.7 counts s 1 in the 20-60 keV energy band 
in the 100-s bin centered around 2003-07-01T20:03:02.552 (in 
pointing 00870035). Meanwhile, the JEM-X flux peaked at 
1 1.4 counts s _1 in the 3-35 keV energy band. 




0.54 

time [UJD-52621] 



Fig.l. ISGRI (I, 20-60 keV) and JEM-X (J, 3-35 keV) 
lightcurves around the flare of IGR J08408-4503 (with a 100 s 
binning), and hardness ratio (I-J)/(I+J). 



The peak 2-10 keV X-ray flux could be estimated as 1.1 x 
10~ 9 erg s _1 cuT 2 , assuming the spectral model described below 
(wabs*comptt). The average flux during the 1000-s flare is 7 x 
10~ 10 erg s _1 cirT 2 . 

Mosaicking 6 pointings surrounding the main flare (3 before 
and 3 after, i.e. excluding the pointing with the flare itself) shows 
that IGR 108408^-503 is also clearly detected around the main 
flare, from 2003-07-01T16:40:00 until 23:35:13. 

A mosaic of all pointing images containing the source (ex- 
cept 8 pointings surrounding the flare) was produced and gave 
a 5-cr upper limit for the flux in quiescence of 0.056 counts s" 1 
in the 23^0 keV energy range, for an effective exposure time 
of 2.45 Ms. The 2-10 keV upper limit can be estimated as 
6 x 10~ 12 



ergs 1 cm 2 . 



The JEM-X and ISGRI average spectra of IGR J08408^1503 
were extracted from the pointing containing the flare, and are 
shown in Fig. [2] The spectrum was fittecj^] with xspec using a 
simple absorbed powerlaw, but no good fit could be obtained 
with xl — 1-7 (for 13 degrees of freedom). Then, two mod- 
els exhibiting a spectral break were used : a broken powerlaw 
(bknpower) and a comptonisation model (comptt). Both mod- 
els give much better results, with the xl for the comptt model 
decreasing to 1.0 (for 12 degrees of freedom). Applying an 
F-test shows that there is a 99.3% probability that the comptt 



2 For all fits, a constant is added to allow for normalisation be- 
tween JEM-X and ISGRI, and a wabs is added to take into account the 
hydrogen column, N H . 



model is better than the simple powerlaw, implying that a spec- 
tral break is needed. The bknpower model suggests a break at 
an energy of 14.8 + 3.4 keV. The comptt model gives a plasma 
temperature kT of 7 + 1 keV, and a plasma optical depth Tpi asma of 
4+1. The intercalibration constant between ISGRI and JEM-X 
is 1.2 + 0.3 in the case of the powerlaw model, and 1.3 ± 0.3 for 
the comptt model. 

Only an upper limit on the value of the hydrogen column Nh 
can be obtained : the 2-cr upper limit on Nh is 3 x 10 23 crrT 2 . 




+ 



chgnnel energy (keV) 



Fig. 2. JEM-X and ISGRI spectra of IGR J08408-4503 (for 
pointing 35 only), fitted with a constant*wabs*comptt model. 



2.2. Swift/XRT observations 



Two Swift ( |Gehrels et al.|2004 1 Target of Opportunity observa- 
tions of IGR J08408^1503 were performed : they consisted of a 
nominal 3.6 and 4.0 ks exposure in photon counting mode start- 
ing at 2006-05-22T10:41:01 and at 2006-07-29T00: 15:00. The 
data were analysed with version 2.3 of the Swift softwar^] 

The best position for IGR J08408^1503 (derived from the 
second observation) is : RA(J2000) = 08 h 40 m 47 s , Dec(J2000) = 
-45°03'32" (with an uncertainty of 4"! at 90% confidenc e 
level), consistent with the results of Kennea & Campana (2006). 

The source is very faint, with only 22 counts in total (resp. 
20 for the second observation), and 18.3 ± 4.7 counts after back- 
ground subtraction (resp. 18.4 ± 4.5). Assuming the spectral 
model derived from INTEGRAL (wabs*comptt), the flux mea- 
sured by Swift can be estimated as 3.4 x 10~ 13 erg s _1 cirT 2 in the 
2-10 keV energy band (resp. 3.1 x 10~ 13 ). However, this model 
is not compatible with the fact that all (but one) events detected 
by Swift have an energy smaller than 1.5 keV. Assuming a soft 
thermal model, one find that these 22 counts are compatible with 
a 0.5-2 keV flux of 9.4 x 10~ 14 erg s" 1 cirT 2 . 



3. Optical counterpart 

3.1. Imaging 

A DSS-II red image of the field is shown in Fig. [3] In the 
2MASS catalogue, there is only one star in the Swift error circle : 
HD 74 194, which has a position of : RA(J2000) = 08 h 40 m 47*79, 
Dec(J2000) = -45°03'30'/2. The next closest source is nearly 
10" away from IGR J08408^503. However, HD 74 194 is so 
bright that we cannot exclude the presence of other faint sources 
in the error circle. Therefore, we cannot conclusively associate 



3 The Swift software is available 
http : //swift . gsf c . nasa . gov/ 



from the GSFC website 
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IGR J08408^1503 with HD 74 194 by positional arguments 
alone. 



HD74194 




54-* 48 42 36 
Right Ascension (J2000) 



30° 



Fig. 3. DSS-II-Red optical image of the field of IGR J08408- 
4503. The outer black and middle blue circles indicate respec- 
tively the INTEGRAL ISGRI and JEM-X error circles, and the 
innermost white circle indicates the Swift XRT position. The star 
HD 74 194 matches perfectly. The field size is 5.8' x 3.7'; North 
is up and East to the left. 



3.2. Spectroscopy 

Spectroscopy of the star HD 74 194 was secured under pro- 
gramme 70.C-(9396(A) on January 9, 2003 at the ESO 3.6m 
telescope equipped with EFOSC, which uses a (2x2 binned) 
2060x2060 pixel Loral/Lesser CCD. Two 5-s spectra were ac- 
quired starting at 08:38:39 UT; grism #12 and a slit width of 1" 
were used, providing a 6 000-10000 A nominal spectral cov- 
erage. The use of this setup guaranteed a final dispersion of 
4.2 A/pix. 

After cosmic-ray rejection, the spectra were reduced, back- 
fire) 



us- 



ground subtracted and optimally extracted (Home 
me Wavelength calibration was performed using He-Ar 

lamps. The spectra were then simply normalized to their con- 
tinuum because no spectroscopic standard was available for the 
star HD 74 194 in the ESO archive, and were stacked together 
to increase the S/N ratio. The wavelength calibration uncertainty 
was ~0.5 A for all cases; this was checked using the positions of 
background night sky lines. 

The analysis of the optical spectrum of HD 74 194 (Fig. |4j> 
shows the Ha line in emission (with EW = 1.0+0.1 A) with a P- 
Cyg profile; besides, Hei absorption lines at 6675 and 7065 A 
and the diffuse interstellar bands (DIBs) at 6 280 A and 6 515 A 
are also detected. 

The P-Cyg profile indicates an escape velo city of 600 + 
50 km s _1 . Assuming that Voo/vgscape = 2.78 ±0.36 (Groenewegen 



et al. 1989 ), this is consistent with the value of Voo = 2000 km s 
derived from IUE spectra by Lamers et al. ( |1995| l and |Snowetal.| 
1994). We also analysed the archival HST/GHRS spectrum of 



HD 74 194 around the O iv and Si iv lines, and found maximum 
velocities of 3100 + 100 kms -1 for both lines (see Fig. [5]). 

The positional consistency of this star with the quiescent 
emission of IGR J08408-4503 suggests that HD 74 194 is an 
O-type supergiant fast X-ray transient. 
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Fig. 4. Left : Combined 6000-9000 A spectrum of the star 
HD 74 194, acquired with the 3.6m ESO telescope. The main 
spectral features are labeled. The symbol © indicates atmo- 
spheric telluric absorption bands. Right : HST/GHRS 1380— 
1410 A spectrum of HD 74 194, showing the S iv and O iv lines 
(datasetZ37J6)2QDT). 



The fact that HD 74 194 is of spectral type 08.5 lb (Wa lborn 
[1973] ) implies My = -6.2 and (B - V) Q = -0.285 ( |Lang|1992| i. 
Usin g the observed V = 7.55 and B - V = +0.22 for this 
star ( |Driliing||1991| l, we obtain a color excess of E(B - V) = 
+0.505 and a distance d ~ 2.7 kpc. This distance estimate 
is in agreement with that obtained from the Hipparcos satel- 
lite (Perryman et al. 1997| l, which measured an annual paral- 
lax 7T = 0.36 + 0.64 milliarcsec for HD 74 194; this value im- 
plies a distance d ~ 2.8 kpc (with a large error) to this source. 
We also estimated the color excess from IUE archival spectra 
(datasets LWP2966[l-2] and SWP53Q2 [3-5]), by calculating 
the reddening necessary to compensate for the 2200 A dust ex- 
tinction feature, and found E(B - V) = +0.63, similar to the value 
derived above when taking into account the uncertainty on the 
extinction law. The color excess is less than half of the value 
of +1.66 obtained from infrared maps (Schlegel et al. 1998[ l, 



thereby supporting a relatively close distance to Earth. This dis- 
tance places HD 74 194 in the Carina Galactic spiral arm tan- 
gent (see, e.g., |Leitch & Vas isht 1998), and in this direction, 
the interstellar absorption column density can be estimated as 
M? M 3 x 10 21 cm" 2 . 



4. Discussion 

Fast hard X-ray outbursts are mostly known in three broad 
classes of objects discussed below : Be transients, thermonuclear 
flashes, and supergiant fast X-ray transient (SFXT) systems. 



4.1. Be transient systems 

The outbursts exhibited by most Be transient systems are of 2 
different origins : they are either related to the orbital period of 
the neutron star reaching the disk of the Be star; or they are large, 
seemingly random, and last a few days. 

Few Be HMXBs are known to show fast and occasional out- 
bursts. An example is A0538-66, which is probably a super- 
Eddington X-ray source ( |Skinner et al.|1982 i. 

The flares of IGR J08408^1503 are probably not of this type, 
because the timescales are completely different, its flaring lumi- 
nosity is less than 1 % of the Eddington luminosity, and its optical 
counter-part is not a Be-type star. 

4.2. Thermonuclear flashes 



4 The image reduction and analysis facility (IRAF) software is avail- 
able from the NOAO website : http : //ir af . noao . edu/ 



Alternatively, thermonuclear flashes - similar to those seen in 
LMXBs - could be envisioned. But in HMXBs, the presence of 
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a magnetic field funnels the accretion stream onto the magnetic 
polar caps of the neutron star (NS), thus preventing the formation 
of a boundary layer in which H and He could burn. Moreover, 
these type I flashes are shorter (=; 30 s), have a much higher peak 
luminosity, and are much softer than observed in IGR J08408- 
4503 (Remilla rd et al.|20 06). Therefore, the origin of the flares 
in this case is most likely not linked to thermonuclear flashes. 

4.3. Supergiant Fast X-ray Transients 

IGR J08408-4503 shares many charact eristics with other S FXT 
systems discovered by INTEGRAL (e.g. Sgu era et al.|2 006) : 



- two X-ray flares have been observed from IGR 108408^-503 
in 12 weeks of cumulated INTEGRAL observations; 

- the bright flares are fast, with a typical duratiorj^Jof 1000 s; 

- the X-ray luminosity during the 1000 s flare is 6 x 
10 35 ergs -1 , while the quiescence luminosity is lower than 
5 x 10 33 ergs -1 and has been as low as 3.0 x 10 32 ergs -1 
during the Swift/XKI observation; 

- during flares, the hard X-ray spectrum features a cutoff at 
14.8 + 3.4 keV that is typical of accreting pulsars; 

- the optical counterpart is a supergiant OB star. 

It is impossible to formally exclude that the X-ray source may 
be unrelated to the O star. However, for the rest of the dis- 
cussion, the many similarities between IGR J08408^-503 and 
other SFXTs will be taken as an evidence for its association with 
HD 74 194. 

The typical X-ray luminosity of O-type stars ranges from 
10 31 to 10 34 ergs -1 , so it is likely that the X-ray luminosity ob- 
served in quiescence (3 x 10 32 ergs -1 ) is emitted either by the 
star itself (through shocks in its wind), or by the interaction of 
its wind with the interstellar medium. The very soft spectrum 
observed in quiescence also does not show any evidence for ac- 
cretion. 

HD 74 194 features micro-variability with a period of 
7.8 days (Koen & Eyer 2002). Such a micro-variability has been 
associated with stellar spots, or with density waves in the atmo- 
sphere ( |van Genderen et al.|1 989). Nevertheless, the variability 
amplitude of 5% does not seem to be related to the observed X- 
ray variability and is too large to be the signature of an orbital 
period. The high luminosity X-ray flares, with their hard spectra, 
are to be explained by another mechanism. 

The stellar wind parameters of HD 74 194, in particular its 
mass loss rate M - 1.4 x 10 -6 Mq yr -1 ( Krticka & Kubat|200i) 
and its terminal velocity v M — 2000 ± 300 kms -1 ( Lamers ~et al.| 
|1995| l, can be used as input for different models. 



If a NS orbits within the wind of an O-type star, an X-ray 
luminosity of 10 36 (r or b/10 12 ) -2 ergs -1 (where r or b is the orbital 
radius in cm) is expected, which is much higher than observed 
on average in the present case, unless r or b is larger than 10 13 cm. 



The flip-flop instability of wind accretion flow (Benensohn 
et al.|1997] > cannot account for the observed variability : it pre- 
dicts variations of the accretion rate by a factor up to 10 only, 
i.e. much smaller than observed, and on a time scale R^/v^ 
10 2 - 10 3 s (where R d is the effective accretion radius), much too 
short to explain the observed flare frequency. 

The propeller mechanism (Illarionov & Sunyaev 1975 ) could 
also be envisioned to explain the outbursts, supposing that during 
quiescence, the NS is rotating fast enough to centrifugally inhibit 
the accretion of matter from the O-type companion. Outbursts 



5 The total duration of the first flare of IGR J08408-4503 observed 
by INTEGRAL was at least 6 hours. 



would then be observed when a change occurs, switching off 
the propeller mechanism for a short period. Considering the ob- 
served burst duration, the changes cannot be driven by rotation. 
It is also very unlikely that the ignition and stop of accretion 
could be driven by a sudden change of the magnetic field and 
back within a few ks. 

The most likely source for the observed flaring behavior is 
therefore a sudden change of the density of the gas trapped 
within the accretion radius : this is expected to happen if the 
stellar wind is clumpy. 

Assuming that the flares are the signature of the accretion of 
wind clumps on a compact object orbiting at 10 13 cm, and that 
the flare duration is related to the clump size, the clump column 
density can be evaluated as 5 x 10 23 cm -2 . Assuming a spheri- 
cal shape for the clumps, their typical mass can be estimated as 
10 18 g. The flare frequency can then be used to estimate the stel- 
lar mass loss rate in the form of clumps as 10 -6 Mq yr -1 (for a 
spherical geometry). 

The density of the clumps (~ 10 13 gem" 3 ) is only slightly 
above the peak densities predicted by hydrodynamic simulations 
of growing wind perturbations induced by line-driven instabili- 
ties (see Fig. 6 in Feld meier et al.|1997 l. The predicted density 
ratio of 10 4 -10 6 between the dense and hot phases of the wind 
at a given radius may explain the large variations between the 
flaring and quiescent luminosities. 

The wind clumps are massive enough to generate the out- 
bursts while being rare enough to explain their low frequency. 
Using the observed luminosities and flare frequencies, the clump 
mass loss rate and Nh density are highly dependent on the or- 
bital radius considered. If this simple model is correct, the flare 
frequency]^] constrains this radius to be of the order of 10 13 cm. 
The main difference between SFXTs and persistent super-giant 
HMXBs would then be their different orbital radius : a few R* in 
persistent systems, but tens of R t in SFXTs. 
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